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Abstract 



It can be found that with the adequate measure, the beta stabihty hue 
C3 , is equidistant from neutron and proton drip lines. We explore this fact 

and its predictive potentiality in the simplest case, the classic liquid drop 
formula. 

^ I For a nucleus of atomic number A = N + Z in the beta stability line, we 

■ can consider the corresponding nuclei {Z — k, N) in the proton drip line and 
(Z, N — k') in the neutron drip line, with respective masses Ap, An. The main 

' mass models in the market (eg, from |2]) predict a very small difference Ap~An, 

, which even becomes zero in isolated points under the action of microscopical 

■ corrections. 

We have studied this difference for the classical Weizsacker formulapp , 

.1/^ {A-2Zf 



u . 

2 ■ An analytical -even if very large- expression can be given if instead of taking 

A as the independent variable, we fix the mass Aq in the drip lines. Then solving 
the second degree equation in the proton drip line 



M[Z, Aq] - M[Z - 1, Ao - 1] - = 

and the third degree one in the stability line (we take nip — m„ ^ but it is 
not necessary) 

we can get the corresponding mass A and proton and neutron numbers Z, N{= 
A — Z) of the stable nucleus. We compare this neutron number with the one 
got from the neutron drip line equation 

M[Ao - No, Ao] - M[A -No,Ao-l]- m„ = 0. 

The difference (i(Ao) = N — Nq results a very convenient function to input in 
a numerical-analytical program such as, for instance, Mathematica, because we 
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can plot dependences with any of the four free parameters of the model, as 
well as mixed plots d{Ao,ai) or d{ai,aj). It is specially relevant to check the 
dependence in 03, because it has a natural minimum for the zero value, but it is 
not uniformly increasing; there is a second minimum in the ^ 1 MeV area, but 
this one has also a dependence on Aq so we can no expect it to coincide exactly 
with the usual value 03 = 0.711 MeV. Still, this minimum can be interpreted as 
the cause of our equidistancy. 
For the usual values 

ai = 15.75MeV,a2 = 17.8MeV,a3 = 0.711MeV,O4 = 23.0MeV, 

it can be seen that d{Ao) gets the higher value for Aq 300; it is only -0.815 
when proton and neutron masses are equal, and this maximum discrepancy 
moves by only about two units when proton and neutron masses are given 
different value, so for simplicity one can find convenient to keep with mp — m„ ~ 
as we do here. 

The function is not linear, so for mid-range masses the difference is appre- 
ciably smaller. Generically we can say that the equidistance property k = k' 
with different proton and neutrons masses holds within a 2%. 

As a numerical experiment, and to test how strong a rule is the equidistance 
to the driplincs, we measured the square discrepancy at four different values, 
using an averaging function 

mean = ^(50)^ -|- ^(100)^ + d{150f + d{200f, 

and we took the bold bet of searching for minima^ in the allowed parameter 
space of ai...a4. 

After discarding unphysical or non- convergent zones, we explored a detailed 
area for 04 = 6. .18, 02 = 0..40, 03 = 0.3. .1.2, 04 = 5. .45. We seed the algorithm 
with values in this area and then we ask it to descend numerically into a local 
minimum. Note that smaller values of 03 cause the algorithm to descent into the 
trivial 03 = case. The results arc seem in the table. The stability line depends 
only of 03/04, so it is a good indication to get both increasing, even if below the 
empirical value. It can be said that we are always out, in the four parameters, 
by a sensible percentage, but if one takes into account that no empirical input 
is used in the calculation^, the exploration is successful enough. Besides, we 
have chosen the averaging function in a very arbitrary way. This choosing is 
surely more important than the seed. For the usual values, d{Af)) has three 
zeros besides the trivial, about A = 25, 600, 2000. This is very generic, and it 
means that up to three discrete values can be adjusted in an unnatural manner. 
The averaging formula must evaluate a number of points enough to be sure that 
no zeros happen in the range of validity of the mass formula, say A = 30... 300. 
Prom the shape of d{Ao) in this range it seems that a weighting proportional 
to the binding energy of Aq (say, in the stability line) could be adequate, but it 

1 Mathematica 5.x users beware: the default method in FindMinimum notices the sum of 
squares and it fails; you must use any of the older methods, 
^except, if you want, the value of the atomic mass unit 
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is still an ad-hoc suggestion, as any function with the same shape would do the 
same role. 

Experiments can be done with variations in only three, two or one param- 
eters, but then the extant ones must be fixed with other methods, beyond the 
purpose of this brief note. Contrary to the others, the parameter 04 seems not 
be very able to vary, as it stops about two or three units from the seed. But 
fixing it to 23.0 docs not improve the match. 

As we have said, the equidistance property can be noticed in most models 
of nuclear masses, and our function d{AQ), or alternatively any measure of the 
discrepancy between k and k', is an interesting parameter to consider when 
studying the properties of a mass model. 

An explanation of this property should be that really some important mass 
dependent effect is enhanced at the drip lines, so this effect forces mass formulae 
to adjust themselves to fit. The effect does not need to happen uniformly in all 
the mass range, it is enough to force the coincidence in isolated points, perhaps 
the ones having strong microscopical corrections, as noted above. 

Perhaps it is related to the extra nuclear stability coming from magicity, 
because an additional phenomena appears when we consider magic numbers at 
the drip linos: we can use nuclear mass to pair Z and N numbers, in such way 
that a nucleus in the proton drip line given a magic (or semimagic) Z number 
correspond to another nucleus of the same mass in the proton neutron drip line, 
having again a magic (or semimagic) N number: 



At neutron dripline, N= 


28-30 


50 


(64) 


82 


f26 


184 


At proton dripline, Z= 


28 


(40) 


50 


(58) 


82 


114 



mean sqr discrepancy 


ai 


02 


as 


a4 


0.00329185 


11 


1382 


18.2911 


0.521987 


18 


9961 


0.00450291 


11 


0895 


19.869 


0.556106 


17 


1047 


0.0173419 


13 


9102 


20.5569 


1.62424 


13 


6503 


0.0158552 


14 


8943 


19.1126 


2.07992 


14 


2993 


0.0144648 


15 


4414 


18.2103 


2.34985 


14 


7491 


0.00267294 


12 


2724 


18.3904 


0.568557 


22 


0129 


0.0021288 


13 


5535 


18.1768 


0.606681 


26 


1328 


0.00175986 


15 


0299 


17.7684 


0.668602 


30 


3896 


0.00149186 


16 


6956 


17.0859 


0.748912 


34 


9357 


0.00108511 


18 


1636 


12.6318 


0.812369 


41 


3469 


0.00110731 


18 


2587 


12.8536 


0.823127 


41 


2093 


0.00113735 


18 


1823 


13.6355 


0.811909 


40 


8982 


0.00212587.. .(*) 


14 


2744 


16.8097 


0.70215 


26 


5332 


0.00214133.. .(*) 


13 


7113 


18.2199 


0.628928 


25.981 



Table 1: Some typical predictions from Mathematica' FindMinimum explo- 
rations of the allowed parameter range for the liquid drop model. To be com- 
pared with the usual ai = 15.75, 02 = 17.8, as = 0.711, 04 = 23.0. Lines marked 
with (*) took 15.75 and 17.8 as seed 
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By the way: the author was too worried about other phenomena and he 
thought the fact here reported were well known. So I must thank M. Asorey by 
asking why the drip lines I kept drawing were always equidistant of the stability, 
and A. Zuker by suggesting to look first in the simplest model, Weizsaecker's 
liquid drop. 
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